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energy resources, which produce power with minimum greenhouse-gas emissions like wind, photo- 
voltaic and micro turbines. Installation of distributed generation at non-optimal location can result in 
various problems such as an increase in system losses and costs, voltage rise and fluctuations, reliability 
and stability problems. Therefore, it is necessary to develop an optimization or heuristic technique 
Keywords: based methodology to identify the optimal placement of distributed renewable generation for a given 
Distributed generation f system that can provide economic, environmental and technical advantages. There are several 
Distributed renewable generation researches that study on the optimal distributed renewable generation location by their imposed 
Power distribution planning : ie : noe - De it 
Optimization methods. constraints and objectives. However, the systematic principle for this issue is still an unsolved problem. 
This paper reviews some of the most popular distributed renewable generation placement methods, 
including 2/3 Rule, Analytical Methods, Optimal Power Flow, Mixed Integer Nonlinear Programming, 
various types of Artificial Intelligent optimization techniques and Hybrid Intelligent System. Each 
methodology has its own features and potential for significantly promote the applicability of 
distributed renewable generation in power systems. 
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1. Introduction 


One of the most important attractions for the studies on the 
integration of distributed resources to the grid is the exploitation 
of the renewable resources, such as wind, solar, hydro, biomass, 
geothermal and ocean energy. The ‘renewable’ term is defined as 
primary, domestic and clean or inexhaustible energy resources [1]. 
They are naturally scattered around the world and also are 
smaller in size. Accordingly, these resources can only be tapped 
through integration into the distribution network by means of 
distributed generation (DG). 

Although there is no consensus on the exact definition of DG, 
there are some significant attempts, in the literature [2,3], to 
define the concept. DG is any source of electric power of limited 
capacity, directly connected to the power system distribution 
network where it is consumed by the end users. Governments 
around the world strive to achieve ambitious targets of incorpor- 
ating considerable amounts of distributed renewable generation 
(DRG) and combined heat and power (CHP) in response to the 
climate-change challenge and the need to enhance fuel diver- 
sity [4]. Fig. 1 shows the annual DRG capacities around the world, 
with estimated increasing from 6000 to 17,000 MW in 2009 to 
2015 [5]. 

A multitude of recent events has created a new environment for 
the electric power infrastructure. They are listed as follows [6]: 


(i) Deregulation of the electric utility industry and the ensuing 
break-up of the vertically integrated utility structure. 

(ii) Public opposition to building new transmission lines on 
environmental grounds. 

(iii) Keen public awareness of the environmental impacts of 
electric power generation. 

(iv) Rapid increases in electric power demand in certain regions 
of the country. 

(v) Significant advances in several generation technologies that 
are much more environmentally benign (wind-electric gen- 
eration, micro turbines, fuel cells, and photovoltaic) than 
conventional coal, oil, and gas-fired plants. 

(vi) Increasing public desire to promote “green” technologies 
based on renewable-energy sources. 
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Fig. 1. Annual renewable distributed generation capacity additions, world mar- 
kets: 2009-2015 [5]. 


(vii) Awareness of the potential of DRG to enhance the security of 
electric power supply, especially to critical loads, by creating 
mini- and micro-grids in the case of emergencies and/or 
terrorist acts, and/or embargoes of energy supplies. 


Panvara et al. revealed that over the period from 1971 to 1995, 
CO, emissions grew at an average rate of 1.7% per year and had a 
faster growth rate, at 2.2% per year for the period towards 
2020 [1]. Table 1 shows an overview of the most important 
emissions related to electricity production based on different 
technologies [7]. Ackermann et al. reviewed five renewable DG 
technologies in the power sector: hydro, wind, photovoltaic, 
geothermal and tidal power. This selection is fairly representative 
of technologies that are important in terms of their potential 
capacity to contribute to a low-carbon world economy. The data 
comprises direct emissions and indirect emissions. Indirect emis- 
sions are emissions that occur during the manufacturing of the 
power unit and the exploration and transport of the energy 
resources. Currently, only hydropower and wind power generate 
significant low-carbon portions of global electricity. Biomass is not 
included in the table, as it is considered CO2 neutral, as the amount 
of COz emitted into the atmosphere when biomass is burned is 
equal to the amount of CO» absorbed during its growth [7,8]. 

The DRG advantages can be best described by dividing them into 
three aspects, which are technical, economic and environmental. The 
technical advantages cover wide ranges of benefit such as efficiency, 
grid reinforcement, power loss reduction, reliability, eliminating or 
deferring the upgrades of the power system, improving load factors 
and voltage profile, thus increased power quality. The economic 
advantages include the reduction of transmission and distribution 
operating cost, reduced health care costs due to improved environ- 
ment, saving the fossil-fuel cost and reduction on electricity tariff. 
Environmental advantages entail the reductions in emission of green- 
house gases (SO2, CO2), reducing sound pollution, and conservation of 
resources for additional use [9,10]. 

Despite its promises, installing the DRG into an electric power grid 
is not a simple plug-and-play problem. Indeed, as well as the 
operation of the DRG itself, it requires a careful consideration for 
the interaction with existing power network with respect to stability, 
reliability, protection coordination, power loss, power quality issues, 
etc. [11,12]. Installation of DRG at non-optimal places can result in an 
increase in system losses, reconfiguration of protection scheme, 
voltage rise and fluctuations, increase in costs, etc. It is necessary to 
develop proper methodologies and tools, which are capable of 
identifying the optimal location of DRG for a given system. These 
methodologies are based on the heuristic techniques or optimization 
programs [13]. 

This paper aims to review the techniques of DRG optimal 
location and sizing in power system environments, which sum- 
marized as shown in Fig. 2. In Section 2, the conventional 
algorithm of DRG placement methods is presented, which include 
the 2/3 Rule [14], Analytical Methods [15-20], Optimal Power 
Flow [21-29], and Mixed Integer Nonlinear Programming 
[30-33]. Section 3 reviews the Artificial Intelligent methods such 
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Table 1 


Comparison of energy amortisation time and emissions of various distributed generation technologies [7]. 
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Technology 


Energy pay back time in months 
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Fig. 2. Optimization placement techniques for DRG. 
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as the Evolutionary Algorithm [34-54], Simulated Annealing [55-57], 
Differential Evolution [58-63], Particle Swarm Optimization [64-70], 
Fuzzy System [71,72] and Tabu Search [73,74]. Section 4 presents 
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possible promising techniques for future use in DRG placement such as 
the Ant Colony Search Algorithm [75,76], Artificial Bee Colony Algo- 
rithm [77-80], Cuckoo Search [81,82], Firefly Algorithm [83,84], 
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Imperialist Competition Algorithm [85-88], and various types of Hy- 
brid System [89-98]. The conclusions are summarized in Section 5. 


2. Conventional optimization techniques for distributed 
renewable generation placement 


This section reviews the conventional optimal placement 
methods of DRG in power systems. These approaches vary 
according to their objectives such as to minimize power losses, 
minimize cost of systems, improve the voltage profile, increased 
the load ability and local marginal price. The mathematical 
algorithm, including the objective function, to solve the 
optimization-based DRG placement problem is discussed in this 
section. 


2.1. 2/3 Rule 


This rule is often used in proper placement of shunt capacitors 
in distribution systems. The 2/3 Rule proposes that the best 
capacitor size is 2/3 of the load which is located at 2/3 of the 
distance from the feeder. It also can be extended to “2/ 
(2N+1)Rule” for N capacitors. For instance, the optimal locations 
for two units with approximately 2/5 capacity for each of them, 
might be located at the 2/5 and 4/5 length of the line [14]. 

This method is suitable for a radial feeder with uniformly 
distributed load, where it is suggested to install DRG of approxi- 
mately 2/3 capacity of incoming generation at approximately 2/3 
of the length of the line. It was proposed for minimizing the losses 
and voltage impacts and for feeders with uniform loads. Since it is 
a simple approximation, this technique cannot be applied directly 
to feeders or to distribution network systems with different types 
of loads such as increasing loads, centrally loads and other real 
systems with non-uniform loads. 


2.2. Analytical Method 


Various Analytical Methods (AM) have been proposed for the 
placement of DRG with optimal size in the distribution network. 
Most of the methods are based on theoretical, mathematical 
analysis and calculations [15-20]. They have common goals, 
which are to reduce the power losses, improve voltage profiles, 
finding optimal size and optimal location. Fig. 3. depicts the AM 
structure in solving the DRG planning problem. 

In a work published in 2004, Wang and Nehrir [18] presented 
AM to determine the optimal location in radial as well as network 
systems to minimize the power loss of the system. To find the 
optimal location of DRG, the objective function (f;) for DRG at bus 
j is as follows: 


j-1 N 
f= XO Rui@|Sul’+ D> Ru |Sul, j = 2,3,- -N (G9) 


i=1 i=1 


where, Rj((j) is the equivalent resistance between bus 1 and bus i 
when DRG is located at busj,j + 1. Sz; is complex power: 
f Real(Z11 +Zä—2Z1i), i<j 
M+ Real (Zig Benen ee i>j 2) 
eal (Z11 +Za-1i-1) 1-9) i>j 


where Z411, Zi, Zi; are the elements of impedance matrix. 
When the DRG is located at bus 1 (j=1), the objective function 
will be as follows. 


N 
fi = DOR |Suil? (3) 


i=1 


The goal is to find the optimal bus m where the objective 
function reaches its minimum value as follows: 


fm =Min(f;), j=1,2,..N (4) 


The theoretical procedure to find the optimal bus to place DRG 
in a network system can be summarized as follows: 


(i) Admittance matrix is calculated without DRG, then admit- 
tance matrix, impedance matrix, and equivalent resistances 
are calculated for different DRG locations. 

(ii) Objective function values for DRG are calculated at different 
buses to find the optimal bus m. 

(iii) If all the voltages were in an acceptable range when the DRG 
is located at bus m, then bus m is the optimal site. 

(iv) If some bus voltages do not meet the voltage rule, then move 
the DRG around bus m to satisfy the voltage rule. 

(v) If there is no bus that can satisfy the voltage regulation 
rule, then try a different size of DRG and repeat the 
procedure. 


Acharya et al. [15] also proposed an analytical expression to 
calculate the DRG optimal size and placement for minimizing the 
total power losses in primary distribution systems. The objective 
function considers the total power loss (P,), as shown in Eq. (5), 
where this expression is popularly known as “exact loss formula” 
in power systems [99]: 


N N 
P= [æy (PiP;+QiQ;) + By (PiP;—Q:Q;)] (5) 
i=1j=1 

where oj = rj/vivj cos(ôi—ô;), By =Try/vivj sin(ði—ð;), rj+Xij=Zij 
are the ijth element of bus impedance matrix (Zpus), Vi, Vj are 
the voltages at ith and jth buses, respectively, P; and P; are the 
active power injection at the ith and jth bus, respectively, Q; and 
Q; are the reactive power injection at the ith and jth bus, 
respectively, N is number of buses, ô; and 6; are voltage phase 
angle at ith and jth buses, respectively. 

Acharya et al. used the concept that approximates loss and 
actual loss calculated by utilizing the accurate load flow, follow- 
ing the same pattern. Using this concept, the load flow analysis 
required only two cycles, first time is by applying to the base case, 
and another time by applying with the DRG attached to obtain the 
final solution. The optimum size of DRG for each bus is calculated 
using Eq. (6), which is obtained by equating the rate of change of 
losses with respect to injected active power from DRG to bus ith 
to zero: 


N 
Poci = Poit+ z BiQi— XO ayPj—bjQ) (6) 
n j=1j+i 

where Ppg; is the real power injection from DRG placed at node i, 
and Pp; is the load demand at node i. Then, the approximate loss is 
computed for each bus by placing DRG with optimal size, one per 
time. The bus corresponding to minimum total loss will be the 
optimum location. After that the load flow analysis with DRG 
gives the final result; the exact loss. 

In 2009, Gözel and Hocaoglu [16] formulated the problem 
using the expression for the total real loss (P,,,;), to find the 
optimum location of DRG and was expressed as a function of the 
branch current injection: 


nb 
Pigss = > |Bi|? -Ri = [R]" |[BIBC] -[1)|? (7) 


i=1 


630 
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Fig. 3. Analytical Method flowchart [15-20]. 


where R; is the ith branch resistance, [R] is branch resistance row 
vector, nb is number of branches, BIBC is bus-injection to branch- 
current matrix, and [J] is the vector of the equivalent current 
injection for each bus except the reference bus. 

The method is based on the equivalent current injection that 
uses the bus injection to branch current (BIBC) and branch current 
to bus voltage (BCBV) matrices, which were developed based on 
the topological structure of the distribution systems and is widely 
implemented for the load flow analysis of the distribution 
systems [100]. The proposed method requires only one base case 
load flow to determine the optimum size and location of DRG. The 
formula is derived as the derivative of the total power losses per 
each bus injected real powers, and equated to zero to determine 
the optimum size. The optimum size DRG is placed at each bus, 
and loss is calculated. The bus corresponding to minimum power 
loss will be defined as the optimum location, if the approximated 
bus voltages are within a limit; otherwise the DRG will be omitted 
from that bus, and next higher loss bus is chosen, and voltages are 
checked for acceptable limit again until the optimum location 
is found. 


In 2010, Hung et al. developed a comprehensive analytical 
expression based on an improvement to the method proposed by 
Acharya et al. [15] in 2006, which was limited to only one DRG 
type, extended to find the optimum sizes and optimal location of 
various types of DRG. Authors considered four major types of DRG 
based, which include DRG capable of delivering real power only, 
DRG capable of delivering both real and reactive power, DRG 
capable of delivering real power and absorbing reactive power, 
and DRG capable of delivering reactive power only [17]. Hosseini 
and Kazemzadeh [20] and Hung and Mithulanantha [19] also 
proposed an improved analytical formulation with the considera- 
tion of the power factor for single DRG placement in order to 
minimize power losses. 


2.3. Optimal Power Flow 


In 2005, Harrison et al. implemented Optimal Power Flow 
(OPF) considering ‘reverse load-ability’ approach to maximize the 
capacity of DRG and identify any available headroom in the 
system within the imposed thermal and voltage constraints. The 
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modelling technique commonly uses the steady-state DRGs as 
negative load [25]. The objective function is as follows: 


n 


fa) = X -Ci x MW? (1-4) (8) 


i=1 


where yw is capacity adjustment factor, MW? is initial active 
power capacity of DRGs in pu, C is capacity value in per unit 
megawatt of DRG capacity, i is DRG bus index, and n is number of 
buses available for capacity addition. 

Vovos et al. dealt with generating the capacity allocation, 
considering additional fault level constraints imposed by protec- 
tion equipment such as switchgear. An iterative process would 
allocate new capacity using OPF mechanisms and readjust the 
capacity to bring fault currents within the specifications of 
switchgear [28,29]. 

Gautam and Mithulananthan [24] formulated the problem of 
optimal placement, including size, with two distinct objective 
functions, namely social welfare maximization and profit max- 
imization. Social welfare is defined as the difference between the 
total benefit to consumers, minus the total cost of production 
[101]. The objective function associated with social welfare has 
been formulated as the quadratic benefit curve submitted by the 
buyer (DISCO), (B(Pp;) minus quadratic bid curve supplied by 
seller (GENCO), C(Ppi) minus the quadratic cost function supplied 
by DRG owner C(Ppci): 


N 
f= X Bi(Ppo)—CiPp)—CPoai)) (9) 


i=1 


The profit maximization is formulated as follows: 


profit; = 2; x (Poci)—C(Ppci) (10) 


where Pp,; is the DRG size at node i; A; is the locational marginal price 
(LMP) at node i after placing DRG; C(Ppci)=dpcit bociPpci)+ 
CpciPpci)* is the cost characteristic of DRG at node i. In this method, 
the traditional OPF algorithm for cost minimization is modified to 
incorporate the demand bids, in addition to the generation bids. The 
candidate locations for DRG placement are identified based on LMP, 
which is determined as the Lagrangian Multiplier of the power 
balance equation in OPF. The base case OPF based on social welfare 
maximizing algorithm is used to evaluate the generation dispatch, 
demand and prices for each of the nodes. The nodal prices obtained 
are used as indicators to identify the candidate nodes for DRG 
placement. The placement is intended to meet the demand at a 
lower price by provide a scenario of variety of DRGs. 

In a work published in 2009, Jabr et al. presented an ordinal 
optimization (OO) method, which aims for loss minimization and 
DRG capacity maximization in siting and sizing of DRG. The 
confidence in the OO solution to solve complex problems is set 
mainly by two parameters; « that specifies the percentage of top 
designs and P that is the probability level, which at least one of 
the sampled designs belongs to the top & percentage. The OO 
approach can be easily implemented by using the existing OPF 
software [102]. Algarni and Bhattacharya [21] integrated the 
goodness factors of DRG units directly into the distribution 
system operation model, based on the OPF framework for incre- 
mental contribution of DRG unit to active and reactive power 
losses, termed as incremental loss indices. 

Dent et al. [23] assessed the maximization of total generation 
under network security constraints using an OPF model which 
was solved by gradually adding limited numbers of line outage 
contingencies, until a solution to the complete problem is 
obtained. The limit on the number of contingencies added 
reduced the size of the optimization problems encountered. 
Moreover, varying the limit provides a highly efficient means of 


searching for multiple local optimal of the nonlinear optimization 
problem. Apart from the above OPF-based method, OPF also has 
been used by authors in [22] for evaluating the maximum 
capacity of variable DRG. Ochoa et al. proposed a multi period 
AC OPF to determine the optimal accommodation of renewable 
DRG, which aims to minimize the system energy losses. The 
algorithm is also applied to evaluate the maximum capacity of 
new variable DRG which is able to be connected in a distribution 
network when active network management control strategies are 
in place. The active network management schemes embedded 
into the OPF include coordinated voltage control, adaptive power 
factor and energy curtailment [26,27]. 


2.4, Mixed Integer Nonlinear Programming 


Mixed Integer Nonlinear Programming (MINLP) refers to mathe- 
matical programming with continuous and discrete variables and 
nonlinearities in the objective function and constraints. In 2005, 
El-Khattam et al. proposed an integrated model for solving the 
distribution system planning problem by implementing DRG as an 
attractive option in distribution utility’s territories. The total invest- 
ment objective function is based on the supply-chain model formula- 
tion. It aims to minimize the investment and operating costs of DRG, 
payments toward purchasing the required extra power by the DISCO, 
total payments toward compensating for system losses along the 
planning period, as well as alternative costs according to the available 
different scenarios. The DISCO may have the following alternatives to 
serve its demand growth [31]. 


(i) Scenario A: Purchasing the required extra power from the 
main grid and pumping it to its distribution network through 
its junction substation with main grid. 

(ii) Scenario B: Purchasing the extra power from an existing inter 
tie and delivering it to its distribution network territory. 

(iii) Investing in DRG as an alternative for solving the distribution 
system planning (DSP) problem without the need for feeder 
upgrading. 


The objective function is as follows [31]: 


N T M 
J= X Cri (SB& +BK) opci+8760 X> >> P'CriSoci 
i=1 t=1i=1 
T Av? 
8760 MS —> -pf-C. 
£8700 3 DMD PC: 
= = j=1 Y 
+ Cost of Scenario A or+Cost of Scenario B (11) 


where 
Cost of Scenario-A is as follows: 


SS TU TN M T 
Ca= » 5 CiuOiu K 2 > Cyoij +8760 5 Bf CeSiu Cin (12) 


i=1lu=1 i=1j=1 t=1 
Cost of Scenario-B is as follows: 
TN M T 


TU 
Ca= X X Cjoyj +8760 X f" > PÍ Cint(Sint) Sint Fine (Sint) (13) 


i=1j=1 t=1 i=1 
and 
Bi =10 +a) (14) 


In the above formulation the factors such as backup DRG unit 
capacity(BK), discount rate (d), investment cost (C), operating 
cost (C,), electricity market price (Ce), cost of feeder (Cj), cost of 
transformer (C;,,), intertie power cost (Cint), number of load buses 
(M), power generated from DRG (Spc), power imported by intertie 
(Sint), transformer u in substation i dispatch power (S;,,), number 
of substation (SS), incremental time interval (t), horizon planning 
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year (T), total number of buses (TN), total number of substation 
transformers (TU), feeder segment impedance (|Zj|), system 
power factor (pf), DRG binary decision variable (opg), feeder i to 
j binary decision variable (oj) transformer u in substation i binary 
decision variable (ciu), intertie binary decision variable (cint), DRG 


capacity limit (sper) are considered. 


The proposed model integrates comprehensive optimization 
model and planner’s experience to achieve optimal sizing and 
siting of distributed generation. In this model, binary decision 
variables are employed in the General Algebraic Modelling 
System (GAMS) [103] where unity or zero decision variables 
would mean to invest or not to invest, respectively, to provide 
accurate planning decisions. 

Subsequently, in 2010, Kumar and Gao [32] used MINLP to 
determine the optimal location and number of DRGs in the pool 
as well as the hybrid electricity market. The main contribution in 
this work is: (i) to find the optimal location of DRG based on real 
power nodal price and real power loss sensitivity index as an 
economic and operational criterion; (ii) to determine optimal 
location and optimal number of distributed generators in the 
identified zone based on the Mixed Integer Nonlinear Program- 
ming based approach; and (iii) to find the impact of demand 
variation for nodal price variation, fuel cost and real power loss 
without and with DRG. The optimization problem has been 
formulated in GAMS using SNOPT solver [103]. MATLAB and 
GAMS interfacings have been used to solve load flow at base case 
to obtain load flow data and other parameters required for 
modelling algebraic equation in GAMS [104]. 

Atwa et al. proposed a probabilistic-based planning technique 
that combines all possible operating conditions of the renewable 
DRG units (i.e. wind-based DRG, solar DRG, and biomass DRG) 
with their probabilities. The problem was formulated as MINLP, 
considering the uncertainty associated with renewable DRG 
sources as well as the hourly variations in the load profile. The 
planning problem comes with an objective function for determin- 
ing the optimal fuel mix of various types of renewable DRG units 
in order to minimize the annual energy losses without violating 
the system constraints. The constraints include the voltage limits, 
the feeders’ capacity, the maximum penetration limit, and the 
discrete size of the available DRG units [30]. 

A two-stage MINLP based methodology was proposed by Pokar 
et al. for DRG placement. The method comes with an objective 
function to minimize the cost and maximize the total system 
benefits (TSB). Optimal placement and size are obtained from the 
total cost minimization mathematical problem, which is solved in 
the first stage. It shows that there are an optimal location and size 
for each DRG cost characteristic and for each investment payback 
time. Consequently, the second stage solves the optimal DRG 
investment payback time results from the TSB maximization 
problem. The various DRG technologies include renewable 
resources like wind and photovoltaic to give the system deciders 
some choices, which offer the opportunity of selecting the right 
energy solution at the right location [33]. 


3. Artificial Intelligent optimization techniques for 
distributed renewable generation placement 


Artificial Intelligent (AI) is commonly defined as the science 
and engineering of making intelligent machines, especially intel- 
ligent computer programs [8]. This section covers a wide range of 
AI techniques such as Evolutionary Algorithm [34-54], Simulated 
Annealing [55-57], Differential Evolution [58-63], Particle Swarm 
Optimization [64-70], Fuzzy Systems [71,72] and Tabu Search 
[73,74], which have been applied in most optimization problems 


as well as DRG optimal placement. The applications and goals of 
these techniques varieties of techniques owe to their great 
potentials to optimize technical and economical DRG challenges. 
The fundamental to these techniques are considered below. 


3.1. Evolutionary Algorithm 


An Evolutionary Algorithm (EA) is meta-heuristic population- 
based optimization process and converge to the global optimum 
solution with probability of one by a finite number of evolution 
steps performed on a finite set of possible solutions [105,106]. 
EA is a subset of evolutionary computation, which includes 
Evolutionary Programming (EP), Evolutionary Strategy (ES), and 
Genetic Algorithm (GA) for optimization based on natural selection: 
crossover, mutation, recombination, reproduction and selection 
operators on the population of individuals to perform the search. 

Crossover involves choosing a random position in the two 
strings and swapping the bits that occur after this position; 
mutation randomly perturbs a candidate solution; recombination 
randomly mixes their parts to form a novel solution; reproduction 
replicates the most successful solutions found within a popula- 
tion; whereas selection purges poor solutions from a population. 
The population is randomly created at the start of the search. 
Fitness is used to select individuals from the current generation to 
advance into the next generation. These individuals are recom- 
bined and possibly mutated to form the next generation. This 
process is continued until there is no change in the best individual 
in the population [107]. 

The EP is introduced first, and followed by ES and GA [105]. 
The feasibility to solve efficiently the optimal siting and sizing of 
DRGs by implementing GA was demonstrated in [45]. GA is a 
searching and optimization method based on a model of evolu- 
tion adaptation in nature. It is a very powerful search algorithm 
and is different from other conventional search algorithms. GA 
does not need derivatives or other auxiliary knowledge. GA works 
with a population of individuals and each individual stands for a 
solution. The quality of a solution is evaluated from its fitness, 
which is calculated by using fitness function. In case of DRG 
placement, this fitness is evaluated based on minimizing real 
power losses, to reduce investments and operational costs, and 
providing optimal size [108]. 

An improved Hereford Ranch Algorithm (HRA) was implemen- 
ted in 1998 with single objective function to minimize the active 
power loss, and compared with second-order method and classi- 
cal GA. From the various test results, HRA indicates superior 
performance over the conventional ones [40]. An EP based 
methodology has been developed and presented by Khatod 
et al. for finding the optimal locations of photovoltaic arrays 
and wind turbine generators in a radial distribution system [51]. 
The active energy loss was minimized, considering the constraints 
on bus voltages, line loadings, number of DRGs to be placed and 
dispatched wind power. Various authors in [42,44,46,48,49] 
implemented GA to handle single objective, which minimize the 
total real power losses in the distribution network. 

Teng et al. [50], Borges et al. [35] and Shaaban et al. [54] 
proposed a strategic DRG placement method by implementing GA 
[35,50,54]. The fitness evaluation function that drives the GA to 
the solution is derived as a benefit/cost relation, where the benefit 
is measured from the reduction of electrical losses, and the cost is 
dependent on investment and installation. GA has also been used 
by Borges et al. to evaluate the DRG impact on reliability, losses 
and voltage profile along with DRG planning [35]. Moreover, in 
2012, Shaaban et al. [54] took into consideration about the 
uncertainty and variability associated with the output power of 
renewable DRG as well as the load variability. 
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GA has been used by many authors to handle multi-objective 
(MO) in DRG placement [34,37,38,47,53]. In 2005, Celli et al. 
proposed a MO formulation for the siting and sizing of DRG 
resources into existing distribution networks. The methodology 
objective was achieved by minimizing different functions, which 
is expressed as 


Min C(X(U)) = Min[C,,C,,Cens, Ce] (15) 


where X(U) is a power flow solution calculated as function of vector 
U, which stores the data about the location and the size of generator. 
Ca is cost of network upgrading, C; is cost of energy loss, Cens is cost of 
energy not supplied, C; is cost of purchased energy. The implemented 
technique is based on GA and an é-constrained method that allows 
obtaining a set of non-inferior solutions [37]. 

Carpinelli et al. [36] extended the MO approach in order to 
include uncertainties in DRG energy production. Each possible 
future is formulated as a scenario [109]. The mathematical 
formulation of objective function is formed with three objectives 
as follows: 


Min C(X(U)) = Min|F,,F2,F3] (16) 


where F; is cost of energy loss, F> is voltage profile, and F; is 
power quality. Subsequently, a “double trade-off method” is used. 
In the first trade-off, by using a MO ¢-constrained technique, 
allows DRG siting and sizing solutions for all the scenarios 
considered; such as being associated to a different set of wind 
speed at all the possible locations; the second one allows isolating 
the most robust solutions. In this way, the planner is completely 
free to drive the optimization in a certain direction without losing 
objectivity and generality [36]. 

In 2007, Haesen et al. proposed a robust MO planning 
methodology for the integration of stochastic generators in 
distribution grids. This approach utilizes the Strength Pareto 
Evolutionary Algorithm (SPEA) algorithm, a predecessor of the 
SPEA2 algorithm [39]. It was updated to the SPEA2 algorithm in 
2009 by Rodriguez et al., extended to include the analysis of 
controllable DER units and enhanced to take account of objectives 
that reflect environmental impact and voltage quality [34]. 
Another MO programming approach based on Non-dominated 
Sorting Genetic Algorithm (NSGA) was applied by Ochoa et al. to 
find the best configuration that maximizes the integration of 
distributed wind power generation while satisfying voltage and 
thermal limit [43]. Kumar et al. [41] in 2010 proposed the DRG 
integration approach with MO model implemented for quick 
restoration and the reduction of additional power demand under 
cold load pickup using GA. In 2011, Moeini-Aghtaie et al. [52] 
implemented NSGA-II to minimize the total costs, total losses, and 
improve system reliability in the distribution system. 


3.2. Simulated Annealing 


Simulated Annealing (SA) is a process in which the optimiza- 
tion problem is simulated through an annealing process. It has the 
ability of escaping local minima by incorporating a probability 
function in accepting or rejecting new solutions. SA was intro- 
duced by Kirkpatrick, Gelatt, and Vecchi in 1983 [110]. Due to its 
implementation together with simplicity and good results, its use 
has been growing since mid-80s [111]. The initial temperature 
and cooling procedure are of paramount importance for the good 
use of SA. The algorithm is based on initialization, perturbation, 
cooling schedule, and acceptance probability to perform the 
search. 

Sutthibun and Bhasaputra [56] presented a model to deter- 
mine the optimal location and size of DRG in order to minimize 
the real power loss (P;), emission (Epg), and the contingency in 
terms of severity index (SI) while being subjected to power 


balance constraint and power generation limit using SA as an 
optimization tool. The MO function (F) is the weighted sum of 
individual objective, expressed as follows: 


F =Wy1P,+W2Epg + W351 (17) 


where Wj, w2, and w3 are weight factors whose values are 
between 0.2 and 0.6 with condition w;+w2+w3=1. 

Another MO optimization model for use to determine the 
optimal sizing and locating DRG was implemented by Aly et al. 
[57] for cost minimization, which was achieved via the mini- 
mization of system losses, complex power acquired from DRG, 
and the number of DRG connected. In 2012, Ghadimi and Ghadimi 
[55] proposed SA in order to minimize power losses for sizing and 
sitting of DRG and capacitor banks in the distribution network. 


3.3. Differential Evolution 


Differential Evolution (DE) developed by Storn and Price [112] 
is a simple yet powerful population based, stochastic function 
minimizer and has been found efficient and effective to solve 
various natures of engineering problems [113-116]. DE algorithm 
is a population-based algorithm like GA, using similar operators; 
crossover, mutation and selection. DE is used to solve the 
optimization problem by sampling the objective function at 
multiple randomly chosen initial points. Pre-set parameter 
bounds define the region from which ‘M’ vectors in this initial 
population are chosen. DE generates new solution points in ‘D’ 
dimensional space that are perturbations of existing points, by 
applying three fundamental operators explained below. 


i. Mutation: DE mutates and recombines the population to 
produce a population of ‘M’ trial vectors. Differential mutation 
adds a scaled, randomly sampled, the vector difference to a 
third vector as follows: 


VP = Xie +0 (X-X) (18) 


base 


where ø is known as scale factor and usually lies in the range 
[0, 1], x” and x are two randomly selected vectors (p + q); 


x is known as base vector; V\“is a mutant vector; the base 
vector index ‘base’ may be determined in various ways. This 
may be a randomly chosen vector (base 4 p # q). 

ii. Crossover: DE employs a uniform crossover strategy. The 


crossover generates trial vectors, i, as follows. 


a if (rand; < C, or j =jrana) 

= $ Vo : (19) 

Xij otherwise 

where C, is crossover probability which lies in the range [0, 1], 

defined by the user which controls the number of parameter 

values copied from the mutant. rand; is the jth evaluation of a 

uniform random number generator with outcome [0, 1]. jana 

is randomly chosen index, which ensures o get at least one 

parameter from v\"). ' 

iii. Selection: Objective function is evaluated for target vector 
F , and trial vector f(X®); trial vector is selected if it 
provides better value of the function than the target vector as 


follows: 


xen i a EEP JE) 20) 


x otherwise 


In a work publish in 2010, Hejazi et al. [63] presented DE for 
optimal allocation of distributed generation in distribution net- 
works considering various technical and economic aspects of the 
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problem and formulate multi-objective function, which include 
the cost of network upgrading, cost of purchased energy, cost of 
energy losses, total voltage deviation and total capacity release. 
Estabragh and Mohammadian [62] also presented DE for finding 
the best location of DRG to reach the best condition of voltage 
stability security margin and power losses. Load ability limit is 
formulated as the voltage stability index for security assessment 
in nonlinear optimization problem. 

Various authors [58-60] also proposed DE for the placement of 
single DRG unit in electrical distribution systems to reduce the 
power losses. Abbagana et al. [59] considered two DRG units with 
three different types of generation in the DRG planning metho- 
dology. In 2012, Arya et al. [61] proposed a technique to quantify 
the technical benefits with use of sensitivity approaches to decide 
the appropriate location of wind turbine DRG with induction 
generator, and determination of their optimum capacities by 
using DE to minimize transmission losses. Incremental voltage 
sensitivity aspect was also considered, which accounted for 
voltage stability consideration. 


3.4. Particle Swarm Optimization 


Particle Swarm Optimization (PSO) is a population based 
stochastic optimization technique developed by Dr. Eberhart 
and Dr. Kennedy in 1995, inspired by social behaviour of bird 
flocking or fish schooling [117]. The PSO, as an optimization tool, 
provides a population-based search procedure in which indivi- 
duals called particles change their position (state) with time. In a 
PSO system, particles fly around in a multidimensional search 
space. During flight, each particle adjusts its position according to 
its own experience (this value is called Lpest), and according to the 
experience of a neighbouring particle (this value is called Gpest), 
made use of the best position encountered by itself and its 
neighbour [108]. 

From mathematical point of view, the above procedures for 
each agent can be stated as following: 


vit! — evi +C *rand* (Loest-X') + Coxrandx (Grese—X') (21) 
xit! =—Vvitl, xi (22) 
The weighting function is usually used as follows: 
Mmax—@ min 
W = Wmax ean Itre (23) 


min and @max are the minimum and maximum weights, respec- 
tively. Appropriate values for @ jn and @max are 0.4 and 0.9. Itre 
and Itremax are the current and maximum iteration, respectively. 
Appropriate value ranges for C4 and C, are [1,2], but 2 is the most 
appropriate in many cases. 

A multi-objective formulation for optimal siting and sizing of 
DRG resources in distribution systems was proposed by Hajizadeh 
and Hajizadeh [64] based on PSO and weight method in order to 
minimize the cost of power losses (C,) and energy not supplied 
(Cens). The fitness function was derived from the objective 
function as follows: 


Min C = W,C,+W2Cegns (24) 
where, 
Wi+W2=1 (25) 


Lalitha et al. [65] proposed a two-stage methodology for the 
optimal DRG placement. In the first stage, fuzzy approach was 
used to find the optimal DRG locations and in the second stage, 
PSO was used to find the size of the DRGs corresponding to 
maximum loss reduction. A hybrid objective function was used 
for the optimal DRG placement by Mohammadi and Nasab [66]. 


It had two parts of objective function; in first part, the power loss 
was proposed as an index named Power Loss Reduction Index 
(PLRI) as consideration; in the second part, the effect of DRG on 
reliability improvement on the system was considered, named as 
Reliability Improvement Index (RII). 

In 2011, El-Zonkoly et al. considered a wide range of technical 
issues, which are active and reactive power losses of the system, 
the voltage profile, the line loading and the MVA intake by the 
grid in DRG planning [70]. The PSO technique was implemented 
by Kansal et al. [69]and Zareiegovar et al. [67] to find the optimal 
size and optimum location for the placement of DRG in the radial 
distribution networks by reduction in real power losses and 
enhancement in voltage profile. Voltage stability was taken into 
consideration by Zareiegovar et al. [67] as an objective function in 
DRG planning. 

In order to evaluate the influence of DRG on the total real power 
losses in the power system and to contribute to the DRG impact 
evaluation vis-a-vis the network fault protection strategies, in 
2012, Maciel et al. proposed a multi-objective based PSO (MEPSO) 
for minimizing the real power loss (ILp) and short circuit level 
(ISC3), which are represented through two indices as follows [68]: 


; LossP¢ 

Min ILp = Lossi (26) 
5 IscabcP¢ 

Min ISC3 = ; max Ge (27) 


where Loss”° is the real power for the network with a given DRG 
configuration; Loss? is the total real power loss without DRG; 
IscabcP© represents the three-phase fault current magnitude in 
node i for the network with a given DRG configuration; Iscabc? 
stands for the three-phase fault current magnitude in node i for the 
network without DRG; NN is the number of nodes. 


3.5. Fuzzy System 


The concept of Fuzzy System (FZ) was introduced by Zadeh 
[118] as a formal tool for dealing with uncertainty and soft 
modelling. The fuzzy set theory is widely being applied in power 
systems [119]. A fuzzy variable is modelled by utilizing a 
membership function which assigns a degree of membership to 
a set. Usually, this degree of membership varies from zero to one. 

The data and parameters used in DRG placement are usually 
derived from various sources with a wide variance in accuracy. 
For example, load is considered as known and specified in almost 
all methods, in spite of having a high uncertainty. In addition, 
electricity market price, cost of DRG, peak power saving, etc. may 
be subjected to uncertainty to some degree. Therefore, uncertain- 
ties due to insufficient information may generate uncertain region 
of decisions. Consequently, the validity of the results from 
average values cannot represent the uncertainty level. To account 
for the uncertainties in the information and goals related to 
multiple and usually conflicting objectives in DRG placement, 
the use of fuzzy set theory may play a substantial role in decision- 
making [107]. 

In 2006, the multi-objective allocation of resources was 
implemented by Ekel et al. using Bellman-Zadeh approach as 
the decision-making in a fuzzy environment. Authors developed 
corresponding Adaptive Interactive Decision-Making System 
(AIDMS) base on the Bellman-Zadeh approach [120]. Lalitha 
et al. and Kumar et al. implemented fuzzy set theory to determine 
the suitable locations for DRG placement. Two objectives were 
considered while designing a fuzzy logic for identifying the 
optimal DRG locations, which were to minimize the real power 
loss and to maintain voltage within permissible limits. Voltage 
stability and power loss indices of the distribution system were 
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modelled into voltage stability index (VSI) and power loss index 
(PLI) to obtain DRG suitability index (DSI) as output [71,72]. 


3.6. Tabu Search 


Tabu Search (TS) is a heuristic algorithm for guiding the search 
to find a good solution for a combinatorial problem. It is derived 
from the works of Glover and Hansen both in 1986 for solving 
combinatorial optimization problems [121]. It is an efficient 
combinatorial method that can achieve an optimal or suboptimal 
solution within a reasonably short time. It does not need many 
iteration counts to obtain better solution. It can eliminate local 
minima to search area beyond the local minima. It is based on 
moves, neighbourhood, tabu list, aspiration, intensification, and 
diversification. TS has been successfully applied to obtain optimal 
or sub optimal solutions to problems such as a timetable, a 
travelling salesman and so on [108]. 

To apply TS algorithm for solving optimal operation of dis- 
tribution network, an initial population, which must meet con- 
straints set, is selected randomly. Tabu list is created through a 
selection of a number of members from an initial population that 
has minimum objective function. New population is created 
based on the mutation and recombination rules. Individuals will 
be ranked in a descending order according to their fitness 
function. A number of individuals that have the best fitness 
function are selected into the next population, and the tabu list 
will be updated and checked for convergence [108]. 

The TS application for finding the optimal allocation of DRGs 
from a view point of loss minimization had been illustrated by 
Nara et al. To simplify the algorithm, the determination algorithm 
of the allocation of DRGs and the search algorithm of the sizes of 
DRGs were disconnected, and a decomposition/ coordination 
technique was introduced to implement TS into the problem 
[74]. In 2006, Golsan et al. implemented the TS to solve DRG 
planning problem, including simultaneous DRG, reactive sources 
and network configuration planning. By solving this problem, the 
installation locations, sizes and operation of DRG resources and 
reactive-power sources in selected buses of a distribution system, 
can be determined, along with tap positions of voltage regulators 
and network configuration. In the algorithm, various memory 
structures such as short, intermediate and long-term memories 
are also been implemented. In this work, forbidden moves are 
introduced into the tabu lists by recording numbers that corre- 
spond exclusively to each forbidden move [73]. 


4. Possible promising techniques for future use in distributed 
renewable generation placement 


In recent years, there have been numerous meta-heuristic 
algorithms implemented to solve the DRG planning problem, 
which consist of Ant Colony System Algorithm (ACSA) [75,76], 
Artificial Bee Colony Algorithm [77-80], Cuckoo Search [81,82], 
Firefly Algorithm [83,84], Imperialist Competition Algorithm 
[85-88] and several kinds of Hybrid System [89-98]. 


4.1. Ant Colony System Algorithm 


One evolutionary method that has recently been considered is 
the implementation of finding shortest path by ants. Ant Colony 
algorithms are based on the behaviour of social insects with an 
exceptional ability to find the shortest paths from the nest to the 
food sources using a chemical substance called pheromone [122]. 
The pheromone is the chemical material deposited by the ants, 
which serves as critical communication media among ants. Ant 
Colony System Algorithm (ACSA) is the extension of the Ant 


Colony Optimization (ACO), and it has a better performance than 
ACO in most engineering applications [123-126]. To select the 
next path for any ant, a state transition probability is defined as 
follows. 


(28) 


where qty is intensity of pheromone between nodes i and j; ny is 
the factor indicating the relation of the objective function of any 
problem with the probability of that path; « and f are control 
parameters for determining the weight of trail intensity and 
length of path. 

In a work published in 2008, Wang et al. developed an 
optimization procedure based on the ACSA to seek out the 
optimal re-closer and DRG locations by minimizing a composite 
reliability index. The reliability is enhanced, and authors sug- 
gested that the idea can be extended to the simultaneous 
placement of both re-closers and DRGs, which are dependent on 
one another [75]. Sookananta et al. [76] proposed ACSA for the 
determination of the optimal location and sizing of DRG in radial 
distribution systems to reduce the overall line losses of the 
network. 


4.2. Artificial Bee Colony Algorithm 


Artificial Bee Colony Algorithm (ABCA) is one of the most 
recently defined algorithms by Dervis Karaboga [127] in 2005, 
motivated by the intelligent behaviour of honey bees. ABCA uses 
only common control parameters such as colony size and max- 
imum cycle number. ABCA was initially proposed for numerical 
optimization [127], and it is then applied for combinatorial [128], 
unconstrained and constrained optimization problems [129,130]. 
ABCA employs only three control parameters (population size, 
maximum cycle number and limit) that are to be predetermined 
by the user. ABCA is quite simple, flexible and robust [131-133]. 

In the ABCA system, artificial bees fly around in a multi- 
dimensional search space and some (employed and onlooker 
bees) choose food sources depending on the experience of 
themselves and their nest mates, and adjust their positions. In 
order to produce a candidate food position from the old one in 
memory, the employed bees use the formula as follows [134]: 


Vmi = Xmi + Pmi(Xmi—Xki) (29) 


where x, is a randomly selected food source, i is a randomly 
chosen parameter index and Pmi is a random number within the 
range [—a,a]. After producing the new food source, Vm, its fitness 
is calculated and a greedy selection is applied between Vm and Xm. 
The fitness value of the solution, fitm(xm), may be calculated for 
minimization problems using the following formula [135]: 


a oe i Xm) > 0 
fitm(Xm) = l 1+fm(Xm) iF Fm m (30) 
1+ abs (fn(Xm)) if fm(Xm) < 0 

where fm(Xm) is the objective function value of solution Xm. 

An artificial onlooker bee chooses a food source depending on 
the probability value associated with that food source. The value 
of Pm is calculated as follows [135]: 

fitm(%m) 

=N a (31) 
e aR 1 fitm(Xm) 

Scouts fly and choose the food sources randomly without using 
experience. If the nectar amount of a new source is higher than 
that of the previous one in their memory, they memorize the new 
position and forget the previous one. Thus, ABCA system com- 
bines local search methods, carried out by employed and onlooker 
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bees, with global search methods, managed by onlookers and 
scouts, attempting to balance the exploration and exploitation 
process [135]. 

Latitha et al. and Hussain et al. proposed ABCA to solve the DRG 
placement problem. The objective function was to minimize the total 
system real power loss subjected to equality and inequality con- 
straints. The validity of the proposed method has been proven from 
the comparison ABCA with other existing methods, such as AM and 
PSO. The results showed that ABCA exhibited excellent solution 
quality, fast convergence characteristics and the potential for solving 
complex power system problems [78,79]. ABCA was also implemen- 
ted by Sohi et al. for loss reduction and line capacity improvement in 
DRG planning problem [80]. In 2011, Abu-Mouti and El-Hawary [77] 
proposed ABCA to solve the mixed integer nonlinear optimization 
problem, to determine the optimal DRG location, size, and power 
factor with the aim of minimize the total system real power. 


4.3. Cuckoo Search Algorithm 


Cuckoo Search Algorithm (CSA) is an optimization algorithm 
inspired by the brood parasitism of cuckoo species, which lay 
their eggs in the nests of other host birds. CSA is proposed by 
Yang and Deb in 2009 [136], and it has been applied into the 
engineering optimization problems and shown its promising 
efficiency. For ease in describing CSA, the three idealized rules 
are described as follows [137]. 


(i) At one time, each cuckoo only lays one egg, and leaves it ina 
randomly chosen nest. 

(ii) The algorithm will carry over the best nests with high- 
quality eggs (solutions) to the next generations. 

(iii) A host bird can discover a foreign egg with a probability, 
Pa € [0,1] while the number of available host nests is fixed. In 
this case, the host bird can either abandon its nest and build 
a completely new nest elsewhere or simply throw the 
eggs away. 


When generating new solutions x(t+1) for a cuckoo i, a Lévy 
flight is performed: 


xt) — x 4 o @ Levy(/) (32) 


where « > 0 is the step size which should be related to the scales 
of the problem of interests; in most cases, & can be set to value 1. 
The above equation is essentially the stochastic equation for 
random walk, which is a Markov chain whose next location or 
status only depends on the current location or status, the first 
term, as in Eq. (32), and the transition is probability, the second 
term. The product ® represents the entry wise multiplication, 
which is similar to those used in PSO. In terms of exploring the 
search space, random walk via Lévy flight is more efficient as its 
step length is much longer in the long run [137]. 

The random step length of Lévy flight, which fundamentally 
provides a random walk, is derived from a Lévy distribution with 
an infinite variance and infinite mean. 


Levy ~u = t~} (33) 


here, the sequential jumps of a cuckoo fundamentally form a 
random walk process with a power law step length distribution 
with a heavy tail. Several new solutions should be generated by 
the Lévy walkabout in finding the best solution, as this process 
will speed up the local search. However, a significant fraction of 
the new solutions should be generated through far field rando- 
mization, and whose locations should be far enough from the 
current best solution, to make sure the system not trapped in a 
local optimum [136]. 


In a work presented in 2012, Fard et al. implemented CSA to 
improve voltage profile and reduce power losses in a prolonged 
period for two types of the most economical types DRG units, 
which consist of biomass and solar-thermal [82]. Intelligent 
building of the distribution network is modelled via non- 
constant loads connected to the network with the aid of Monte 
Carlo method. CSA was also applied by Moravej et al. for optimal 
DRG allocation to reduce power loss and improve voltage profile 
of the distribution network [81]. The voltage profile which is the 
main criterion for power quality improvement is expressed by 
two indices, which are voltage deviations from the target value, 
which must be minimized, and voltage variations from the base 
case of network without DRG, which must be maximized. 


4.4. Firefly Algorithm 


The Firefly Algorithm (FA) is a meta-heuristic algorithm, 
inspired by the idealized behavior of the flashing characteristics 
of fireflies. The primary purpose of a firefly’s flash is to act as a 
signal system to attract other fireflies. The idea is to use the 
concept of flashing characteristics of fireflies to consequently 
develop firefly-inspired algorithms as the following three rules 
[138]: 


(i) All fireflies are unisex so that one firefly is attracted to other 
fireflies regardless of their sex. 

(ii) Attractiveness is proportional to their brightness, thus for 
any two flashing fireflies, the less bright one will move 
towards the brighter one. The attractiveness and brightness 
will both decrease as their distance increases. If no one is 
brighter than a particular firefly, it moves randomly. 

(iii) The brightness or light intensity of a firefly is affected or 
determined by the landscape of the objective function to be 
optimized. 


In other words, the brightness can be defined in a similar way 
as the fitness function in GA [138]. In the FA, there are two 
important issues of the variation of light intensity and the 
formulation of the attractiveness. For simplicity, it is assumed 
that the attractiveness of a firefly is determined from its bright- 
ness which in turn is associated with the encoded objective 
function of the optimization problems. In the FA, there are three 
important issues: [139] 

Attractiveness: The main form of attractiveness function 
f(r)can be by any monotonically decreasing functions such as 
the following generalized form: 


Br) =B,e™ (m>1) (34) 


where r is the distance between two fireflies, f(r) is the attrac- 
tiveness at r=O and y is a fixed light absorption coefficient. 

Distance: The distance between any two fireflies i and j at x; 
and x; is the Cartesian distance as follows: 


d 
ri = xix | = | XO (xix—xik)? (35) 
k=l 
where x;, is the kth component of the ith firefly, x;, x; is the kth 
component of the j-th firefly, xj. 
Movement: The movement of a firefly, i is attracted to another 
more attractive (brighter) firefly j, is determined using 


xitla xi 4 BoM (xixi) +a(rand—0.5) (36) 


where the second term is due to the attraction, while the third 
term is a randomization with « being the randomization para- 
meter and rand is a random number generator uniformly dis- 
tributed in range of 0 to 1. 
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In 2012, Saravanamutthukumaran et al. applied FA to find the 
optimal allocation of DRG to minimize real and reactive power 
losses, line loading, short circuit level index, Mega Volt Ampere 
(MVA) intake by grid, and improve voltage profile for different 
load models [84]. The load models included in the work consisted 
of constant load, residential load, commercial load, and industrial 
load. Sulaiman et al. [83] also presented an application of FA in 
determining the optimal location and size of DRG in distribution 
power networks with the aim to minimize real power losses. 


4.5. Imperialist Competitive Algorithm 


Imperialist Competitive Algorithm (ICA) was proposed by 
Atashpaz and Lucas [140]. This method is a new socio- 
politically motivated global search strategy that has recently been 
introduced for dealing with different optimization tasks. Same as 
the Evolutionary Algorithms, ICA commences with an initial 
population of P countries, which are generated randomly within 
the feasible space. The best countries in the initial population are 
selected as the imperialists and other counties are known as the 
colonies of these imperialists. To build initial empires, colonies 
are divided among imperialists based on each imperialist’s power, 
which is define as normalized cost as follows: 


C" = c"—max (c) (37) 


where cp is the cost of n-th imperialist and Cn is the normalized 
cost. As a result, the normalized power of each imperialist can be 
defined as follows: 


(38) 


After dividing the colonies among the imperialists, these 
colonies will start being closed to their empire. The total power 
of an empire can be determined by expression as follows: 


Power, = cost(imperialists,)—é*meancost(colonies of empiresn) 
(39) 


where Power, is the total power (cost) of the nth empire and é is a 
positive number less than 1 and near to 0. Competition among 
empires will then occur, which results as an increase as the power 
of strong empires and decrease as the power of weak ones. Weak 
empires will lose their power and collapse. Finally, by using the 
movement of colonies towards their relevant imperialist and also 
the collapse mechanism, one imperial will exist, which is the 
optimum solution. 

Jahani et al. [88] applied ICA to minimize the total real power 
losses in the distribution network. ICA was also used by Nejad 
et al. [87] to minimize the total real power loss and improve 
voltage profile. In order to maximize the benefits of distribution 
network operators (DNO), Soroudi and Ehsan [85] proposed ICA in 
DRG planning problem and considered the sum of active loss 
reduction and network investment deferral incentives as the 
objective function to be maximized. The optimal location and 
size of DRG units within the network are found considering 
various techno-economic issues, and the objective function is 
formulated as follows: 


OF = p1— Hn (40) 
where 
Np 
[4 =Wx (ies KEDD rr) (41) 
i=1 


Np 
mais XPE (42) 
i=1 
w is the coefficient of incentive for each MWh reduction of active 
losses, Loss” Sis the active loss when no DRG unit is installed in 
the network, Pir is the net active power injected to bus i, y is the 
coefficient of incentive for each MW of installed DRG units, and 
pes is the active power injected by a DRG in bus i. 

In 2012, Rahmatian et al. [86] proposed ICA to determine the 
capacity and location of DRG by considering the islanding mode of 
a distribution network including sensitive loads. The distribution 
system is divided into zones, and each zone should include a DRG 
unit and sensitive loads. On the other hand, when the load density 
is high, the zone will be less extensive and enclose fewer numbers 
of buses. Therefore, the capacity of the installed DRG units will be 
in an acceptable range for each zone. As a result, sensitive loads 
are fed more reliably into the islanding mode, and the network 
losses are reduced. 


4.6. Hybrid Intelligent System 


Hybrid Intelligent System (HIS) is defined as a combination of 
methods and techniques from artificial intelligence sub fields, in 
parallel, to form a new software system. Various types of HIS 
which include Genetic-Fuzzy (GA-FZ) [90,93], Genetic-Tabu (GA- 
TS) [89], Genetic-Particle Swarm (GA-PSO) [94], Genetic-Optimal 
Power Flow (GA-OPF) [91,92,98], Particle Swarm Optimization- 
Optimal Power Flow (PSO-OPF) [96], and Tabu-Fuzzy (TS-FZ) [95] 
have been applied in DRG allocation and sizing problems. 


4.6.1. Genetic-Fuzzy 

Kim et al. presented Genetic-Fuzzy (GA-FZ) method to solve 
DRG placement for distribution systems. Power loss cost of 
distribution systems was taken as the objective function, and 
the number or size of DRGs and deviation of bus voltage were 
taken as constraints. The original objective function and con- 
straints were transformed into equivalent multi-objectives func- 
tions and modelled with fuzzy sets to evaluate their imprecise 
nature. The authors obtained a global solution of multi-objectives 
and imprecise information using goal programming and GA, 
without any transformation for this nonlinear problem to a linear 
model or other methods [93]. 

Haghifam et al. proposed a MO model consisting of monetary 
cost index, operation costs, planning costs, technical risk and 
economic risk. A fuzzy approach was used in the modelling of 
load and electricity price uncertainties and related risks. To solve 
this MO problem, the concept of Pareto optimality, based on Non- 
dominant Sorting Genetic Algorithm (NSGA-II) [141], was used. 
The final solution was found using a max-min approach to select 
the best Pareto-optimal DRG placement solution. The obtained 
results showed that the proposed DRG placement model was a 
powerful decision-making tool for risk management in distribu- 
tion networks with DRG installation and operation [90]. DRG 
planning problems were also solved using GA-FZ by Akorede et al. 
[97] to maximize the system loading margin as well as the profit 
of the DISCO over the planning period. The objective functions are 
fuzzified into a single multi-objective function, and subsequently 
solved using genetic algorithm (GA). 


4.6.2. Genetic-Tabu Search 

A new implementation of a Genetic-Tabu Search (GA-TS) 
algorithm into the optimal allocation of DRGs in the distribution 
network had been illustrated by Gandomkar et al. The power 
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losses, including harmonic power losses were considered as the 
objective function. The proposed algorithm effectively solved the 
DRG allocation problem as demonstrated through the study 
example. A comparison was made between GA-TS, and the GA. 
GA-TS was found to be much better than GA in terms of solution 
accuracy and convergence process. The theoretical procedure of 
the proposed algorithm can be summarized as follows [89]. 


(i) The variables of GA and TS are initialized. 

(ii) An initial population is created by randomly generating a set 
of feasible solutions (chromosomes). 

(iii) Each chromosome is evaluated by running the load flow 
program to determine the fitness function of each chromo- 
some in the population. 

(iv) GA operators are applied to generate new populations as 
follow: copy the best solution from the current into the new 
population, using TS algorithm, generate new members in 
the new population (typically 3-15% of the population size) 
as neighbours to randomly selected solutions in the current 
population. 

(v) The crossover and mutation operator are applied to produce 
new populations. 

(vi) If the convergence criterion is satisfied, stop. Otherwise, 
repeat the process from the step of evaluate the fitness of 
each chromosome. 


4.6.3. Genetic-Particle Swarm Optimization 

In 2011, Moradi et al. present a novel Genetic-Particle Swarm 
Optimization (GA-PSO) for determining the optimal location and 
sizing of DRG on distribution systems. The objective was to 
minimize network power losses, better voltage regulation and 
improve the voltage stability within the frame-work of system 
operation and security constraints in radial distribution systems. 
The proposed methodology was a searching technique developed 
for optimal sitting and sizing of DRG. The problem consisted of 
two parts. The first was the optimal location of DRG, which was 
solved using GA and the second was the optimal sizing solved in 
PSO. GA was chosen to play this role because of its attractive 
quality. The answer obtained from GA solution was used in PSO 
algorithm to optimize the sizing for DRG. The PSO had fast 
convergence ability which is a great attractive property for a 
large iterative and time consuming problem. The intera- 
ction between the two algorithms can be summarized as follows 
[94]: 


(i) Time counter is set to 0 and chromosomes, which represent 
that initial candidates sitting of DRG are generated 
randomly. 

(ii) Each chromosome and optimal sizing of DRG is evaluated 
using PSO. 

(iii) Particle population is initialized, matrix and the size of DRG 
is modified. 

e The objective values which are the total real power losses, 
and the voltage profile improvements are modified. 

e Objective function is recorded as the best candidate among 
the particle and the minimum value as the current overall 
global best of the group. 

e The velocity (v) and position are updated. 

e The stop criterion is checked; if it is satisfied, then stop. 

v) The time counter is updated, t=t+1. 

v) New population of sitting of DRG is created by repeating GA 
(selection, crossover and mutation) until the new population 
is completed. 

(vi) Fitness is evaluated using PSO, and time is updated. 


(i 


(vii) The stop criterion is checked, if it is satisfied, then stop, or 
else go for time updating. 


4.6.4. Genetic-Optimal Power Flow 

Harrison et al. had emphasized that GA combined with OPF 
(GA-OPF) could provide the best combination of sites within a 
distribution network for connecting a predefined number of 
DRGs. In doing so, it overcomes known limitations inherent in 
current available techniques to optimize the DRG capacity. Its use 
would be to enable DNOs to search the network for the best sites 
to strategically connect a small number of DRGs among a large 
number of potential combinations [91,92]. 

For capacity evaluation using OPF, the linear function is 
assumed as follows: 


forr = > Ce (Pe) —Ci (PEM —PA ) (43) 
g=1 


where C,(P,) is the benefit or incentive (£EkW~! year-1) of con- 
necting a generator g of capacity Pz, C, represents the value of the 
loss incentive as applied to the difference between the actual 
level of losses P<’ and the target losses P“. The GA generates 
combinations of locations from those available in the network. 
For each combination of locations, an OPF is performed to define 
the available capacity and evaluate the objective function. This 
information is fed back to the GA which searches for the ‘DNO 
optimal’ connectable capacity. In doing this the method should 
deliver the best locations as well as the capacities available for a 
user-specified number of DRG. 

An OPF was proposed by Naderi et al. in 2012 to minimize the 
capital costs for network upgrading, operation and maintenance 
costs, and the cost of losses for handling the load growth for the 
DRG planning horizon [98]. The term “dynamic” is used to refer to 
the planning over a specific period so that the dynamic distribu- 
tion system planning is, in fact, proposed, where a year- 
dependent decision-making variable would be defined into the 
model. The load duration curve was used to include the way that 
the customer loads may change. Also, the impact on the elec- 
tricity market was considered based on a load-dependent elec- 
tricity price. Moreover, a modified GA was used to find the 
optimal topology solution. 


4.6.5. Particle Swarm Optimization-Optimal Power Flow 

A new hybrid method, which employs discrete PSO and OPF, is 
proposed by Gomez-Gonzalez et al. to connect within a distribu- 
tion network a predefined number of DRG among a large number 
of potential combinations. In PSO-OPF, a particle is defined as a 
combination of DRG locations. The objective function includes 
loss minimization and DRG capacity maximization at once. In 
each iteration and for each particle, an OPF is performed, which 
generates the optimal capacity for DRG unit. This complies with 
the technical constraints imposed by DISCO and agrees to a value 
of the objective function. Computer simulation has demonstrated 
good performance, accuracy, robustness, diversification and 
intensification of the proposed method in comparison to GA [96]. 


4.6.6. Tabu-Fuzzy 

Ignacio et al. present a new possibility’s model by implement- 
ing FZ for the MO optimal planning of power distribution net- 
works that seek for the non-dominated MO solutions 
corresponding to the simultaneous optimization of the fuzzy 
economic cost, level of fuzzy reliability, and exposure (optimiza- 
tion of robustness) of such networks, using an original and 
powerful meta-heuristic algorithm based on TS. This Tabu-Fuzzy 
(TS-FZ) model is used to determine the optimal location and size 
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of the future feeders and substations in the distribution networks 
with dimensions significantly larger than in other papers. The 
model also allows determining the optimal reserve feeders which 
include the location and size that provide the best distribution 
network reliability at the lowest cost for a given level of robust- 
ness [95]. 


5. Conclusions 


This paper provides a timely review of the state-of-the-art of 
optimal DRG placement techniques. Each method has been 
utilized to solve problems with various and limited objectives 
and constraints. Fig. 4 shows a summary of citation, a summary of 


E no. of citations 


20 

18 Piee 

16 

14 

12 

10 

8 

6 

4 

SST kiii 

0 a a a rT m a 
S 2 $ $ K S S $ s R $ P S S KA S 


Fig. 4. Summary of citation of researches on DRG planning in this article. 
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potential researches related to DRG planning, clearly indicating 
the increasing trend over recent years. The techniques imple- 
mented in this literature of DRG allocation are summarized in 
Tables 2 and 3, leading to the following conclusions: 


(i) The analytical techniques have the best precision factor while 
power flow and Artificial Intelligent methods have next 
advance levels. Moreover, the process speed in the Artificial 
Intelligent techniques is slower, compared to the Optimum 
Power Flow and Analytical Methods. However, the considera- 
tion of nonlinear algorithms and integer variables will make 
the conventional method be possibly less robust and the 
running time would be much longer, while Artificial Intelli- 
gent methods can address the integer variable very well. 

(ii) Genetic Algorithm is competent of getting a solution near 
global minima computationally intensive. The Ant Colony 
System Algorithm, Artificial Bee Colony, Cuckoo Search and 
other future promising algorithms still have not been paid 
much attention by researchers. Fuzzy set theory is suitable to 
model the uncertainties in objective function, constraints, 
electricity price, generation, and load for better compromised 
solution. Tabu Search is an efficient combinatorial method 
which can be used to achieve an optimal or sub optimal 
solution within a reasonably short duration. 

(iii) Artificial Intelligent methods are more heuristic than con- 
ventional techniques, which needs further improvement and 
investigation regarding performance on various larger power 
systems with their improved versions, such as the Hybrid 
Intelligent Systems. 


In practice, the main problem is the complexity of the power 
system. Many constraints should be considered simultaneously 
such as power loss, reliability, load factors, voltage profile, voltage 


Brief comparison on the advantages and disadvantages of main approaches applied for the DRG placement. 


DRG placement Advantages Disadvantages Reference 
method 
2/3 Rule Simple and easy to use approximate technique Not suitable applied directly to non-uniform load network [14] 
systems 
Analytical Method Easy to implement, high precision factor, computational time Fewer literature examples, lacks robustness, only can consider [15-20] 
efficiency single objective and single DRG at a time 
Optimal Power Easy to find literature examples, high precision factor, computational Problem formulated in “closed” manner, hard to include [21-29] 
Flow time efficiency different aspect into calculation 
Mixed Integer High precision factor, computational time efficiency Hard to implement and understand [30-33] 
Nonlinear 
Programming 
Evolutionary Efficient performance for finding the global optimum, easy to find Relatively harder to code, premature convergence, possibility [34-54] 
Algorithm literature examples of trapping into local optima, lower precision factor 
Simulated Ease of implementation, ability to provide reasonably good solutions Relatively lower performance for finding the global optimum, [55-57] 
Annealing for many combinatorial problems, robustness large computational time 
Differential Fewer parameters setting, capable of handling complex optimization Unstable convergence, possibility of trapping into local optima [58-63] 
Evolution problems 
Particle Swarm Easy to code with few equations, easy to find literature examples Relatively lower performance for finding the global optimum, [64-70] 
Optimization fewer literature examples 
Fuzzy System Easy to understand and suitable to model uncertainties for better Fewer literature examples [71, 72] 
compromised solution 
Tabu Search Efficient performance for achieve an optimal or sub optimal solution, Relatively harder to code due to many parameters to be tuned, [73, 74] 
capable to escape from local minimum lower precision factor 
Ant Colony Search Easy to understand and code Probability distribution changes by iteration, uncertain time [75,76] 
Algorithm to convergence, fewer literature examples 
Artificial Bee Capable of handling complex optimization problems, easy to code Fewer literature examples [77-80] 
Colony Algorithm 
Cuckoo Search Easy to code, less parameters setting Slow convergence, fewer literature examples [81,82] 
Firefly Algorithm Easy to understand and code Slow convergence, fewer literature examples [83,84] 
Imperialist Capable of handling complex optimization problems Relatively harder to code due to many parameters to be tuned, [85-88] 
Competition fewer literature examples 
Algorithm 
Hybrid Intelligent Efficient performance for finding the global optimum, capable of Relatively harder to code, fever literature examples [89-98] 


System 


handling complex optimization problems 
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Table 3 


Summary of DRG allocation methods. 
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Author(s) Year Methods Objective(s) Power system, load level DRG References 
location 
H.L. Willis 2000 2/3 Rule Minimize power losses and voltage impact None Single [14] 
C. Wang et al. 2004 AM Minimize real power losses 6 and 30 bus system, Time- Single [18] 
varying and Time-invariant 
load 
N. Acharya et al. 2005 AM Minimize real power losses 30, 33, and 69 bus system, Single [15] 
One load level 
T. Gözel et al. 2009 AM Minimize real power losses 12 bus system, One load Single [16] 
level 
D.Q. Hung et al. 2010 AM Minimize power losses 16, 33, and 69 bus system, Single [17] 
One load level (Mix 
sources) 
R.K. Hosseini et al. 2011 AM Minimize power losses 33 and 69 bus system, One Single [20] 
load level (consider 
power 
factor) 
D.Q. Hung et al. 2012 AM Minimize power losses 33 bus system, One load Single [19] 
level (consider 
power 
factor) 
G.P. Harrison et al. 2005 OPF Maximize capacity of DRG 9 bus system, Multi load Multiple [25] 
level 
P.N. Vovos et al. 2005 OPF Maximize the new generation capacities and energy export 12 bus system, One load Single [28] 
level 
P.N. Vovos et al. 2005 OPF Maximize the new generation capacities and energy export 12 bus system, One load Multiple [29] 
level 
D. Gautam et al. 2007 OPF Maximize social welfare and profit 9 bus system, One load level Single [24] 
A.A.S. Algarni et al. 2009 OPF For disco owned: Minimize cost of active and reactive power 18 and 69 bus system, One Multiple [21] 
from substation bus, cost of active and reactive power from load level 
DRG, and cost of DRG active and reactive power in conjunction 
with goodness factor 
For investor owned: Minimize cost of active and reactive 
power from substation bus, cost of energy purchased from 
DRG, and cost of DRG active and reactive power in conjunction 
with goodness factor 
CJ. Dent et al. 2010 OPF Maximize total DRG active power capacity 73 bus system of consists of Multiple [23] 
three area which contains 
24 buses each. One load 
level. 
CJ. Dent et al. 2010 OPF Maximize total DRG active power capacity 10 bus system of U.K. Multiple [22] 
distribution System. One 
load level 
L.F. Ochoa et al. 2010 OPF Evaluate the maximum capacity for variable (renewable) 16 and 61 bus system of Single and [26] 
generation under a range of ANM schemes including U.K. generic distribution multiple 
coordinated voltage control, adaptive power factor control, system. Multi load level 
and energy curtailment. 
L.F. Ochoa et al. 2011 OPF Minimize energy losses 61 bus system of U.K. Multiple [27] 
generic distribution system. 
Multi load level 
W. El-Khattam et al. 2005 MINLP Minimize cost of investment and operation of DRGs, losses 9 bus system. One load level Multiple [31] 
cost and cost of purchasing power by DISCO from grid 
A. Kumar et al. 2010 MINLP Minimize total fuel cost of DRG and conventional generators, 24 bus system. One load Single and [32] 
and line losses level multiple 
Y.M. Atwa et al. 2010 MINLP Minimize annual energy losses 42 bus system. Variable Multiple [30] 
load level (Mix 
sources) 
S. Pokar et al. 2010 MINLP Minimize cost and maximize total system benefit 30 bus test system. Multi Multiple [33] 
load level 
J.O. Kim et al. 1998 HRA Minimize real power losses 6, 14, and 30 bus system. Multiple [40] 
One load level 
D.K. Khatod et al. 2012 EP Minimize active energy losses 69 bus system. One load Multiple [51] 
level 
A. Silvestri et al. 1999 GA Minimize sum of cost of power losses, network reinforcement 43 and 93 bus systems. One Single [45] 
and energy production cost load level 
J-H. Teng et al. 2002 GA Maximize the DRG benefit cost ratio 40 bus system. One load Multiple [50] 
level 
N. Mithulananthan et al. 2004 GA Minimizing real power losses 30 bus system. One load Single [42] 
level 
D.H. Popovic et al. 2005 GA Maximize DRG capacity 75 bus system. Multi load Multiple [44] 
level 
C.LT. Borges et al. 2006 GA Maximize the DRG benefit cost ratio 39 bus system. One load Single [35] 
level 
D. Singh et al. 2007 GA Minimize real power losses 16, 37, and 75 bus system. Single [46] 
Multi load level 
R.K. Singh et al. 2009 GA Minimize real power losses Single [48] 


One load level 
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Table 3 (continued ) 
Author(s) Year Methods Objective(s) Power system, load level DRG References 
location 
30 bus system. Multi load 
level 
R.K. Singh et al. 2009 GA Minimize performance indices include real power losses, 16 and 37 bus system. Single [49] 
reactive power losses, line power flow, and node voltage Constant, residential, 
industrial and commercial 
load models 
D. Singh et al. 2009 GA Minimize performance indices include real power losses, 16 and 37 bus system. Single [47] 
reactive power losses, line power flow, and node voltage Constant, residential, 
industrial and commercial 
load models 
A.A.A. El-Ela et al. 2010 GA Minimize multi-objective function includes voltage profile 8 bus system. One load level Single [38] 
improvement, spinning reserve increase, power flow reduction 
in critical lines and total line losses reduction 
V. Kumar et al. 2010 GA Minimize multi-objective function includes load to be 33 bus system. One load Multiple [41] 
curtailed, bus voltage violation, branch current violation, level 
substation transformer over loading 
M.F. Shaaban et al. 2011 GA Minimize cost of energy losses and cost of interruption 38 bus system. Mixed load Multiple [54] 
level 
N.C. Yang et al. 2012 Dual GA Maximize allowable DRG capacity within short circuit Taipower radial network Multiple [53] 
capacity limits 
G. Celli et al. 2005 GA ands- Minimize multi-objective function include cost of network 78 bus system. Peak load Multiple [37] 
constrained upgrading, energy purchased, energy losses and energy not with constant growth rate 
Method supplied 
G. Carpinelli et al. 2005 GA and £- Minimize multi-objective function include cost of energy 18 bus system. Peak load Single [36] 
constrained losses, improvement in voltage quality and harmonic with constant growth rate 
method distortion 
E. Haesen et al. 2007 SPEA Minimize total line losses, main grid energy flow, DRG 30 bus system. One load Multiple [39] 
installation costs and gas distribution grid investments level 
A.A. Alarcon-Rodriguez 2009 SPEA2 Minimize multi-objective function includes annual DRG UKGDS radial network. Single [34] 
et al. dispatched energy for local ancillary, annual DRG curtailed Stochastic load 
energy, CO, emission, and voltage quality index 
L.F. Ochoa et al. 2008 NSGA Multi-objective functions are to maximize energy export, 33 bus system. Time Single [43] 
minimize real power losses, and minimize single phase short varying load (variable) 
circuit level 
M. Moeini-Aghtaie et al. 2011 NSGA-II Minimize total costs, total losses, and improve system 27 bus system. One load Multiple [52] 
reliability level 
T. Sutthibun et al. 2010 SA Minimize multi-objective function includes power loss, 33 bus system. One load Multiple [56] 
emission, and severity index level 
A.I. Aly et al. 2010 SA Minimize power losses, complex power acquired from DRG, 68 bus system. One load Multiple [57] 
and the number of DRG connected leve 
N. Ghadimi et al. 2012 SA Minimize power losses 33 bus system. One load Multiple [55] 
level 
H.A. Hejazi et al. 2010 DE Multi-objective functions include minimize cost of network 37 bus system. One load Multiple [63] 
upgrading, cost of purchased energy, cost of energy losses, leve 
total voltage deviation and total capacity release 
MRR. Estabragh et al. 2011 DE Minimize power losses and maximize load ability limit 14 bus system. One load Multiple [62] 
level 
J. Gunda et al. 2011 DE Minimize power losses Indian Electricity Board 25 Multiple [58] 
bus system. One load level 
I. Hussain et al. 2012 DE Minimize power losses 33 bus system. One load Multiple [60] 
leve 
M. Abbagana et al. 2012 DE Minimize power losses 33 bus system. One load Multiple [59] 
leve 
L.D. Arya et al. 2012 DE Minimize power losses 6 and 30 bus system. One Multiple [61] 
load level 
A. Hajizadeh et al. 2008 PSO Minimize the cost of power losses and energy not supplied 33 bus system. One load Multiple [64] 
eve 
M.P. Lalitha et al. 2010 PSO Minimize real power losses 33 bus system. One load Multiple [65] 
eve 
M. Mohammadi et al. 2011 PSO Minimize real power losses and improve system reliability 12 bus system. One load Single [66] 
level 
S. Kansal et al. 2011 PSO Minimize power losses 33 bus system. One load Single [69] 
level 
A.M. El-Zonkoly 2011 PSO Minimize multi-objectives indices include real power losses, 38 and IEEE 30 bus system. Multiple [70] 
reactive power losses, line loading, short circuit level, and Constant, residential, 
improve voltage profile industrial and commercial 
oad models 
G. Zareiegovar et al. 2012 PSO Minimize power losses, improve voltage stability, and voltage 69 bus system. One load Multiple [67] 
profile level 
R.S. Maciel et al. 2012 MEPSO Minimize real power losses and short circuit level 34 and 123 bus system. One Multiple [68] 
oad level 
M. P. Lalitha et al. 2010 FZ Minimize real power losses 33 bus system. One load Multiple [72] 
evel 
K. Injeti et al. 2011 FZ Minimize real power losses and improve the voltage profile 12, 33, and 69 bus system. Single [71] 
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Author(s) Year Methods Objective(s) Power system, load level DRG References 
location 
K. Nara et al. 2001 TS Minimize real power losses 28 sections and 78 sections. Multiple [74] 
Multi load level 
M.E.H. Golshan et al. 2006 TS Minimize cost function including cost of power losses at peak 33 and 69 bus system. Peak Multiple [73] 
load time, cost of fuel served for energy losses and the cost of load and multi load level 
reactive sources 
L. Wang et al. 2008 ACSA Optimize multi-objective function consists of SAIFI and SAIDI 39 and 394 bus system. One Multiple [75] 
load level 
B. Sookananta et al. 2010 ACSA Minimize real power losses 12 and 15 bus system. One Single [76] 
load level 
M.P. Lalitha et al. 2010 ABCA Minimize real power losses and improve voltage profile 33 bus system. One load Multiple [78] 
level 
F.S. Abu-Mouti et al. 2011 ABCA Minimize real power losses 69 bus system. One load Multiple [77] 
level 
M.F. Sohi et al. 2011 ABCA Minimize real power losses and improve line capacity 33 bus system. One load Multiple [80] 
level 
I. Hussian et al. 2012 ABCA Minimize real power losses 33 bus system. One load Multiple [79] 
level 
M.M. Fard et al. 2012 CSA Minimize power losses and improve voltage profile 33 bus system. Non- Multiple [82] 
constant load (generated by 
Monte Carlo) 
Z. Moravej et al. 2013 CSA Minimize real power losses and improve voltage profile 69 bus system. One load Multiple [81] 
level 
S. 2012 FA Minimize multi-objectives indices include real power losses, 38 bus system. Constant, Multiple [84] 
Saravanamutthuku- reactive power losses, line loading, short circuit level, and residential, industrial and 
maran et al. improve voltage profile commercial load models 
M.H. Sulaiman et al. 2012 FA Minimize real power losses 69 bus system. One load Multiple [83] 
level 
R. Jahani et al. 2011 ICA Minimize real power losses 34 and 69 bus system. One Single [88] 
load level 
H.C. Nejad et al. 2011 ICA Minimize real power losses and improve the voltage profile 33 bus system. One load Multiple [87] 
level 
A. Soroudi et al. 2011 ICA Minimize power losses and maximize network investment 69 bus system. One load Multiple [85] 
deferral incentives level 
M. Rahmatian et al. 2012 ICA Minimize power losses considering the islanding 33 bus system. One load Multiple [86] 
modedistribution network level 
K-H. Kim et al. 2002 GA-FZ Minimize real power losses cost 12 bus system. Multiload Multiple [93] 
level 
M-R. Haghifam et al. 2008 GA-FZ Multi-objective function is in three groups:(i) cost (ii) 9 bus system. Aggregated Multiple [90] 
technical risk (iii) technical risk.Minimize multi-objective load level 
function include: (i) Cost of energy losses, investment cost of 
DRG units; operation and maintenance cost; (ii) substation 
loading, line loading; voltage. (iii) cost of power purchased 
from grid, cost of power generated by DRG 
M.-F. Akorede et al. 2010 GA-FZ Maximize the system loading margin and the profit of the 6 and 30 bus system Single [97] 
DISCO 
M. Gandomkar et al. 2007 GA-TS Minimize power losses 13 and 34 bus system. Single [89] 
Multiple load level 
M.H. Moradi et al. 2011 GA-PSO Minimize network power losses, improve voltage regulation 33 and 69 bus systems. One Single [94] 
and voltage stability load level 
G.P. Harrison et al. 2007 GA-OPF Maximize DRG capacity 69 bus system. One load Multiple [91] 
level 
G.P. Harrison et al. 2008 GA-OPF Maximize incentive to DNO by optimizing DRG capacity and 69 bus system. One load Multiple [92] 
loss reduction level 
E. Naderi et al. 2012 GA-OPF Minimize investment cost, operation and maintenance cost, 9 bus system. Time-varying Multiple [98] 
and cost of power losses load 
M. Gomez-Gonzaleza 2012 PSO-OPF Maximize DRG capacity and loss reduction 30 bus system. One load Multiple [96] 
et al. level 
IJ. Ramirez-Rosado et al. 2004 TS-FZ Minimization of the economic cost, expected nonsupplied 182 bus real system. Multi Multiple [95] 
energy and exposure load level 


stability, operational cost, emission of greenhouse gases, and 
related capacity. These are the main reasons most researchers 
have divided the problem into limited parts along with imposed 
constraints and tried to propose their own solutions. However, 
the systematic principles for DRG placement issue is still an 
unsolved problem despite many optimization methods have been 
proposed. In this literature, several popular methods to select the 
DRG site have been considered along with their potentials and 
applications. Each methodology has its own features and poten- 
tial to significantly promote the applicability of DRG in power 
systems. 
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